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We study the late-time behaviour of a dynamically perturbed rapidly rotating black hole. Consid- 
ering an extreme Kerr black hole, we show that the large number of virtually undamped quasinormal 
modes (that exist for nonzero values of the azimuthal eigenvalue m) combine in such a way that the 
field (as observed at infinity) oscillates with an amplitude that decays as 1/t at late times. This is 
in clear contrast with the standard late time power-law fall-off familiar from studies of non-rotating 
black holes. This long-lived oscillating "tail" will, however, not be present for non-extreme (presum- 
ably more astrophysically relevant) black holes, for which we find that many quasinormal modes 
(individually excited to a very small amplitude) combine to give rise to an exponentially decay- 
ing field. This result could have implications for the detection of gravitational-wave signals from 
rapidly spinning black holes, since the required theoretical templates need to be constructed from 
i ' linear combinations of many modes. Our main results are obtained analytically, but we support 

^ the conclusions with numerical time-evolutions of the Teukolsky equation. These time-evolutions 

provide an interesting insight into the notion that the quasinormal modes can be viewed as waves 
trapped in the spacetime region outside the horizon. They also suggest that a plausible mechanism 
for the behaviour we observe for extreme black holes is the presence of a "superradiance resonance 
cavity" immediately outside the black hole. 
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The generic response of a black hole to dynamical perturbations was first investigated more than thirty years 
^ ' ago. Vishveshwara [Q] was the first to observe, in simple numerical experiments where a Schwarzschild black hole was 
^ pelted by Gaussian wave-pulses, the "ringing" associated with the exponentially damped quasinormal modes (QNMs). 
qh, Shortly thereafter, the late time power-law fall-off (that all perturbative fields decay as t^^'-^^ in the Schwarzschild 
I I geometry) was discovered by Price [Q. Since these first studies, a considerable body of work has established the 
^jpj importance of these two phenomena for black hole physics (see |^ for a review) . 

• • We now know that the late-time tail is due to radiation backscattered by the weak gravitational potential in the far- 

zoneQ. It is independent of the central object's strong-field features (such as the presence of an event horizon), and 
will consequently be the same for (say) a neutron star and a black hole with the same gravitational mass. In contrast, 
• the QNMs represent radiation scattered in the strong-field regime. Hence, they provide the means for making a clear 
I distinction between a neutron star and a black hole. A neutron star has many families of fluid pulsations modes, an 
observation of which would provide useful constraints on the supranuclear equation of state In addition, there is 
a set of modes that arise because gravitational waves can be temporarily trapped by the spacetime curvature caused 
by the presence of the star. These modes share some qualitative features with the QNMs of a black hole, which 
arise as waves are temporarily trapped in the vicinity of the black hole's gravitational potential peak (located near 
the unstable photon orbit at r = 3il/ in the Schwarzschild case). The modes are analogous to scattering resonances 
in quantum physics and thus decay exponentially in time. In a detailed study. Leaver |Q identified the QNMs as 
complex- frequency poles of the relevant Green's function (see also [Q). In addition, he showed that the power-law 
tail, that dominates the black holes dynamical response once the QNMs have died away, originates from the presence 
of a branch cut in the Green's function (customarily placed along the negative imaginary frequency axis). 

Most previous studies of QNMs and power-law tails have focused on non- rotating black holes, and our understanding 
of this case has reached a mature level. The same can not be said about the rotating case, however. A few years ago, 
the QNMs had been calculated also for Kerr black holes ||, ^ , but there were no actual calculations demonstrating the 
presence of power-law tails. Several recent developments have improved our understanding of dynamical rotating black 
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holes. Of particular relevance has been an effort to develop a reliable framework for perturbative time-evolutions of 
Kerr black holes [^0[ |ll, 12, 13|. There has also been recent efforts to analytically approximate the late-time power-law 
tails for Kerr black holes ^ 0. The main conclusions of these studies support the standard picture: The QNM 
signal (typically dominated by the slowest damped mode) gives way at late times to a power-law fall-off. Although 
the QNM frequencies and the late-time power-law are slightly altered by rotation, the results are qualitatively similar 
to the Schwarzschild ones. For example, for a massless scalar field Ori and Barack ^ predict a fall-off (for non- 
extreme Kerr black holes) proportional to t-'-l™l-3-g where q = 0/1 for even/odd I + m. This should be compared 
to the Schwarzschild result t~^^~^ . The difference can be attributed to the rotation-induced coupling of the various 
Z-multipoles in the signal. 



B. Motivation: The observability of QNM signals 

With a new generation of gravitational-wave detectors due to reach their projected sensitivities within the next 
few years, the question whether we can realistically hope to do "black-hole spectroscopy" by detecting QNM signals 
following (say) the formation of a black hole in a supernova is of obvious interest. In principle, such a detection should 
enable us to infer the black-hole mass and angular momemtum and therefore provide an unambiguous identification 
of astrophysical black holes. 

For slowly rotating black holes this presents a serious challenge. To demonstrate this in a simple way we consider 
the gravitational wave signal associated with a certain QNM. Far away from the black hole the associated flux follows 
from the standard formula 

167r ^ ^ 47rr^ dt ^ ^ 

Using this together with 

(2) 

dt 2Te ^ ' 

where Te is the e-folding time of the QNM, and assuming a monochromatic wave (of frequency /) such that h — 2tt fh 
(although we note that this assumption is not justified for rapidly damped QNMs), we get 

eV'^ 1 



h=[-] ^ (3) 

\TeJ TTrf 

Finally, we estimate the effective amplitude achievable after matched filtering as 

Kg = hVN = hy/f7e (4) 

where iV is the number of detected cycles of the signal. Parameterising this result we arrive at (cf. similar estimates 
for pulsating stars §) 



J0 



where S is the radiated energy as a fraction of the black hole's mass M. For a Schwarzschild black hole the frequency 
of the radiation depends on the black-hole mass as / « 12(M0/M)kHz. Given these relations, and recalling the 
estimated sensitivity of the generation of ground based detectors that is under construction (LIGO, VIRGO, GEO600 
and TAMA), we see that it is going to be difficult to see QNM signals from slowly rotating solar- mass black holes. 
The mode-signals are basically too rapidly damped. But there are reasons not to despair. First of all, the situation is 
much more favourable for low-frequency signals from supramassive black holes in galactic nuclei and detection with 
LISA, the space-based interferometric gravitational-wave antenna. Secondly, there have been recent indications that 



"middle weight" black holes, with masses in the range 100 — lOOOM©, exist 17, |18|. For such black holes the most 
important QNMs would radiate at frequencies where the new generation of ground based detectors reach their peak 
sensitivity (~ 100 Hz). If there are such black holes in our galaxy, and they are dynamically perturbed by some 
external agent, we may hope to detect the resultant QNM ringing in the future. 

It has been suggested that QNM signals from rapidly rotating black holes would be easier to detect than what the 
above estimate suggests. This belief is based on the fact that some QNMs become very long lived as a — > M (where 
< a < M is the rotation parameter of the black hole). QNM calculations predict the existence of an infinite set 
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of essentially undamped modes in the extreme Kerr limit ||, |T9j , and previous investigations into the detectability 
of QNM signals have focused on the slow damping of these modes 2^, ^ . It is easy to see how the decreased 
damping of the modes may increase the detectability considerably. But it is also easy to understand that a slowly 
damped mode is only easier to detect than a short-lived one if the modes are excited to a comparable amplitude. This 
then forces us to address (rather difhcult) issues regarding QNM excitation |2^. In particular, we must investigate 
whether it is easier to excite a slowly damped QNM than a short-lived one. Intuitively, one might expect this not 
to be the case. In similar physical situations the build-up of energy in a long-lived resonant mode takes place on a 
time-scale similar to the eventual mode damping. Consequently, it ought to be quite difhcult to excite a QNM that 
has characteristic damping several times longer than the dynamical timescale of the excitation process. This argument 
suggests that the amplitude of each long-lived mode ought to vanish in the limit a ^ M when the e-folding time of 
the mode increases dramatically. Indications that this is the case were provided by Ferrari and Mashhoon some years 
ago H. 

To illustrate this further, we connect (|) to the amplitude of the QNM under consideration. We do this by assuming 
that the mode-signal can be represented by 

h w sin(27r/t) . (6) 

r 

Then we can use to infer the "amplitude" A that corresponds to a certain "total energy radiated through the 
mode" E. We get 



E = ^IfllA^ , (7) 



and the effective gravitational-wave amplitude 



-3/10 A 

h.n^^A^[l-i,) (8) 



Here we have used an empirical approximation for the frequency and damping rate of the fundamental I — m = 2 
mode |2^. When the effective amplitude is expressed in this form it is quite clear that a decrease in A can easily 
compensate for the increase in Te that occurs as the black hole spins up. In other words, in order to correctly discuss 
the detectability of the Kerr QNMs one must necessarily investigate this balance. This provides the prime motivation 
for the present work. 



C. A brief summary of previous work 

In a previous paper we discussed some issues concerning the long-lived QNMs of a Kerr black hole. Our results 
essentially followed by studying a suitably defined excitation coefficient (defined in terms of the asymptotic behaviour 
of the perturbed field, see Section HE of the present paper), A°^*/a„. Provided that one can calculate this quantity 
(which is far from trivial |^) for a given QNM as the black hole spin varies one can estimate the associated change in 
level of excitation of the mode. This follows since we can express the effective amplitude (^) associated with a single 
QNM (with complex frequency aj„) as 



/ Re LUn A°^' 

2 V Im iLJn OLn 

Using results derived in the main part of this paper we find, for an extreme Kerr black hole, that 

where (5 is a positive real constant (see Section IID for its definition) , and the least damped modes correspond to large 
values of the integer n. From this we immediately deduce that the "effective detectability" is exponentially small for 
these modes. 

Similar evidence is provided by recently obtained numerical data |26j relevant for the full range of a. Before we 
discuss these results we recall that the QNMs of a non-rotating black hole are located symmetrically with respect to 
the imaginary frequency axis in the complex tj-plane. For non-zero a this symmetry is broken and co- and counter- 
rotating modes will be affected in different ways. The counter-rotating modes remain rapidly damped, while the 
co-rotating modes become very slowly damped as a ^ M. These conclusions follow from eg. Figure 3 in S. 
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In Figure |l| we compare the "effective gravitational-wave amplitude" obtained from (|9|) for the slowest damped co- 
and counter-rotating QNMs of a Kerr black hole perturbed by a massless scalar field. From the figure we see that, 
even though the co-rotating mode is much longer lived, its "detectability" decreases significantly as a — > M. In fact, 
our results suggest that the effective amplitude of the slowest damped QNM for a rapidly rotating black hole may be 
as much as three orders of magnitude smaller than the amplitude of the corresponding mode of a Schwarzschild black 
hole. This is clearly very bad news as far as detecting the various QNMs of a rapidly spinning black hole is concerned. 

FIG. 1: An assessment of the "detectability" of Kerr black hole QNMs as a — > M. The effective amplitude h^s (as estimated 
from Eq. ^) is plotted as a function of the spin parameter a. We compare the slowest damped co-rotating QNM to the slowest 
damped counter-rotating one. The figure shows that the co-rotating mode (which becomes very long lived in the near extreme 
case) has a much smaller excitation coefficient than the counter-rotating mode as a ^ M. 

However, one further complication must be accounted for. Our argument shows that each individual QNM has a 
small amplitude in the limit of a rapidly spinning black hole. But we need to account for the fact that a large number of 
modes approach the same limiting frequency as a — > AI . Conceivably, these modes could interfere constructively and 
result in a considerable signal. Indeed, our previous estimates suggest that, despite the fact that the individual 
long-lived modes are excited to a negligible level their collective contribution to the late-time signal is significant. This 
then leads to a late-time behaviour of a rapidly rotating black-hole that is not at all well described in the standard 
terms (as a signal mainly consisting of the leading QNM and the late-time power-law tail) . Instead, for an extreme 
black hole the long-lived modes completely dominate the late-time behaviour, by giving rise to an oscillating signal 
with an amplitude that decays as 1/t. In this paper we present the complete derivation of this result, and discuss the 
extent to which it remains relevant also for astrophysical (non-extreme) black holes. 

D. Organisation of this paper 

The remainder of the paper is organised as follows. Part II contains (Section IIA & IIB) the formulation of the 
initial-value problem in the Kerr geometry, a brief description of the analytic properties of the relevant Green's 
function in the complex frequency plane (Section IIC), and a discussion of the long-lived Kerr QNMs (Section IID). 
Part II then concludes with the calculation of the late-time field due to these modes, which leads to the main result of 
this paper (Section HE). Section III concerns numerical time-evolutions and QNM signal reconstructions that amplify 
and strengthen the analytical results. In Part IV we discuss a physical interpretation of the results, and introduce the 
concept of a "superradiance resonance cavity" . Our conclusions are in Part V. Some technical details are discussed 
in Appendices A-C. In Appendix A we present approximate solutions of the Teukolsky equation for near extreme 
black holes first obtained by Teukolsky and Press In Appendix B the analytical properties of the solutions 

of the Teukolsky equation are examined in some detail. Finally, in Appendix C, the realistic case of gravitational 
perturbations is briefly discussed. Throughout the paper geometrised units G = c = 1 have been adopted. 

II. THE CAUCHY PROBLEM IN KERR GEOMETRY 

In the last decade several results have emphasized the surprising accuracy of black-hole perturbation theory in 
situations where one would intuitively expect it to fail. The most celebrated example of this is the close-limit 
approximation for black hole collisions devised by PuUin and Price p^ . It has become clear that perturbative 
methods provide not only a benchmark test for fully nonlinear numerical relativity, but that it often offers a reliable 
alternative or complement [ p9| . In fact, in many problems black- hole perturbation theory captures a significant part 
of the physics, although it obviously will not provide any insights into the non-linear behaviour (unless taken to higher 
orders [pO[). 

Perturbative studies are attractive because of their relative transparency. Once we make the assumption that 
the perturbative field is weak in the sense that its contribution to the spacetime curvature can be neglected, the 
black hole evolution equations can be cast in the form of a wave equation with a complicated effective potential 

|3l|. This makes the initial- value problem amenable to an analytical (or comparatively simple numerical) treatment. 
Furthermore, in studies where the emphasis is on the qualitative behaviour one can make further simplifications by 
considering a massless scalar field. Since the more realistic cases of gravitational or electromagnetic fields are governed 
by perturbation equations very similar to the scalar wave equation, one would expect the simple toy problem (recall 
that no scalar fields have not yet been observed in Nature) to provide generally relevant information. The advantage 
of the scalar field problem is simply that one does not have to go through the several steps involved in recreating 
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(say) the perturbed metric from the employed variables. Consequently, we mainly focus on the dynamics of a massless 
scalar field in the Kerr background in this paper. But since some of our conclusions may prove relevant for future 
gravitational- wave detection we discuss the gravitational field problem in Appendix C. The main conclusions drawn 
from this Appendix is that all the results from the scalar field problem remain relevant also in the gravitational case. 



A. A formal solution 

The standard formalism for studying Kerr black holeperturbations is based on the Newman-Penrose approach. As 
was first shown by Teukolsky almost thirty years ago [p2| , the evolution of any perturbative field can be described by 
a single "master equation" , which in the case of a scalar field takes the form 



(r2+a2)2 



sin^ e 



AMar 9^$ 



df^ A dtdip 
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A sin^ e 



dr 



1 d 



sine 80 v"''' 



(11) 



2Mr + a^. The two 



where we have used standard Boycr-Lindquist coordinates. In (|ll| ) we have also used A = ' 
horizons of the black hole then follow from A = and can readily be found to be r± = M ± V-M^ — 0? . In order to 
simplify (nTl) we introduce the tortoise coordinate , which is defined as 



Integrating, we have for a < M 



and for a = M 



dr 



2Mr+ , / r 
m — 



2A/r_ , / r 
In — 



„, / 7- 
— + 2Af ln( — - 1 



(12) 



(13) 



(14) 



Furthermore, the axisymmetry of the Kerr background allows us to immediately separate out the dependence on the 
aziniuthal angle tp by introducing 



(15) 



The function ^m(t,r,d) then satisfies 
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(16) 



where G — rA/{r^ + a?Y . 

In order to proceed further we use a Laplace-type integral transform; 
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— oo-\-ic 



(17) 
(18) 



where c is some (positive) constant, to bring the problem into the frequency domain. If the initial data for the field 
^m{t,r,9) is given by 



a$„,(i,r, e) 



$n = 



dt 



(19) 
(20) 
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we then get the transformed equation 

Here we have defined the two differential operators 



V 



(r2 + a2)2 r d 



where K = (r^ + 0^)10 — am, and 
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and the source term is 
5($o,$o) 



iuj{r^ + a ) — AiamMr 
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$0 



2\2 



sin^ I 



$0 



(21) 



(22) 



(23) 



(24) 



Finally, we would like to separate the radial and angular dependencies. To do this we assume that we can find a 
function S{9) such that 



CS + ES = Q 



(25) 



where E is an eigenvalue that ensures that 5* is regular at both 6 = and 9 = tt. The angular eigenfunctions are 
then defined via a Sturm-Liouville problem for E. The theory of such problems tells us that the eigenfunctions form 
a complete, orthogonal set on the interval < < tt for each combination of auj and m. Since the functions must 
limit to the standard spin-weighted spherical harmonics (once the factor e*""^ is included) as acu — > it makes sense 
to label each eigenfunction by the integer / (with / e [0, 00]). In other words, we denote the solutions to (|25| ) as 
Sim{i^, 0). We further require that the eigenfunctions are normalised in such a way that 



SmOSlmSnmdO — Sin 



(26) 



with the bar denoting complex conjugation. By assuming that the left- and right-hand sides of (|T]) are both expanded 
in the complete set of angular functions, i.e. that $„i = '^iZ\m\ ^imi^T'i")Sim{0,a,uj) etcetera, we now get 



(PRlr, 



+ (2amw - a'^uj^ - E)A dG ^2 



where 



(r^ + a?Y 



Rlr. 



(r^ + a^)- 



5i™(0,aw)5($o,$o)sinM0 



We note that the effective radial "potential" is 

^, , + (2amu - a^u? - E)A dG 



G' 



(r2 -I- a2)2 dr* 
A solution to ( p7| ) can be obtained via the Green's function technique. Specifically, we find that 

A(r')G«,r, 



i?/m(w, r) 



-Si,n{oj,r')drl 



(27) 



(28) 



(29) 



(30) 



where the primed variables represent the source point and the unprimed represent the location of the observer. The 
required Green's function can be expressed in terms of two linearly independent solutions to the homogeneous version 
of ( p7| ) . These solutions are defined by their behaviour close to the event horizon and at spatial infinity in such a way 
that 



as r — s- r+ , 



(31) 
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and 



Here 



e *'^''* as r ^ rj 



(32) 



^ = ^ ^T77 — = - mtj+ , (33) 



2Afr, 

where a;+ is the angular velocity of the event horizon. Finally, we have 

1 r u"P(r,)u'"(r:) , < < r 



and the scalar field follows after inverting the integral transform: 

*™(i,r,0) = — ^ / i?,™(w,r)5,™(0,aw)e-''^*dw . (35) 

In this way we have arrived at a formal solution to the initial value problem for the Teukolsky equation. 

B. The asymptotic approximation 

Our aim now is to analyse the Kerr problem further and assess the level of excitation of the various QNMs. It is 
generally accepted that this is a far from trivial problem, see [0, for a detailed discussion, and we will not try to 
provide a complete solution here. Instead, we will study a particular model scenario that allows us to make further 
simplifications. Thus we use the so-called "asymptotic approximation" that was first introduced in i.e. assume 
that i) the observer is situated far away from the black hole, and ii) the initial data has support only in the far zone 
(but inside the observer). The first of these assumptions corresponds to >> M, while the second means that 
M << << r*. Given these assumptions we can use the asymptotic behaviour ( ^l| ) and ( ^2[ ) of the solutions to the 
Teukolsky equation to construct the required Green's function analytically. 

The asymptotic approximation is convenient, but is it realistic to introduce these two assumptions? The assumption 
that the observer is located at large distances from the black hole should be relevant for most problems of astrophysical 
interest (and we note that one could readily devise an analogous approximation for the case when the observer is 
located close to the horizon). Furthermore, there is certainly a class of problems for which the second assumption 
holds. Namely, when well-defined wave pulses fall onto the black hole. This is is a standard scenario considered 
in qualitative studies of black-hole dynamics, but it is not a scenario likely to give rise to detectable astrophysical 
gravitational waves. For the astrophysical scenarios that would seem the most relevant, gravitational collapse to form 
a black hole and the coalescence of two black holes following binary inspiral, the assumption that the "initial data" 
is located far from the horizon will not hold. But again, one could devise a similar approximation in which the data 
has support only close to the horizon. Of course, we know from previous studies of related problems [Q that the main 
QNM excitation arises from data in the region near the peak of the perturbative potential (which is located near the 
unstable photon orbit). For such data the asymptotic approximation will not be accurate, and it is therefore likely to 
be of limited use in realistic scenarios. 

Anyway, once we introduce the asymptotic approximation, we get 

piuir, r /lout > 

G«,n)«-^|-^e--+e--| . (36) 
Then, the solution to the radial Teukolsky equation follows from 

Ri^{cj,r) « - / > -n^r^'"''- • 5,™(u;,r')dr: . (37) 



2iuj J^^ [(r')2-Fa2]2 [ A 
and expression (|3^) for the field becomes 

J-oo+^c 2zw 



^out 



duj (38) 
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Here we have defined the two integrals 

/+00 A(r') 
_ J^^J^T^^^'^'^-S.^i-^^'y^ (39) 

Given specific initial data, the calculation of the field requires evaluation of the frequency integral in (^8|) . If we make 
further simplifications by assuming static initial data, i.e. $o ~ Oi we arrive at 

/•TT /' + 00 

I^^ico) ^iLO desmeSimie,Lo) e^'^<W,0)dr', (40) 

Jo J -oo 



C. Working in the complex frequency plane 

In practice, the frequency integrals in (^) and (|3^) cannot be directly evaluated analytically. It is a straightforward 
task to evaluate them numerically for any acceptable initial data, but this does not provide much insight into the 
qualitative features of the solution. As is by now well known, a useful alternative is to analyze the problem in the 
complex w-plane 0|. To do this, we bend the integration contour into the lower half of the complex frequency plane 
and use Cauchy's theorem to calculate the integral as a sum over the contributions from the deformed contour plus a 
sum of the residues of the poles encircled by the contour. Naturally, we need to understand the analytical properties 
of the Green's function in order to follow this approach. 

The properties of the Schwarzschild black hole Green's function have been exhaustively investigated, but the Kerr 
case has not attracted nearly as much attention js^ . We find that we must consider two different cases depending on 
whether we are interested in an extreme (a = M) or a non-extreme (a < M) black hole. The analytical properties 
of the Green's function for the latter case are exactly as in the Schwarzschild problem: An infinite set of simple 
poles, corresponding to the QNMs, are located in the lower half plane and a branch point (originating from the u"p 
function) at w = leads to a branch cut that is usually placed along the negative imaginary axis, cf. Figure 1 of 
. The existence of this branch cut is related to the asymptotic behaviour of the effective potential in the Teukolsky 
equation. Namely the fact that the potential (beyond the centrifugal term) falls off as ln(r*)/rij as +oo, see |3J] 
for a detailed discussion of this point. 

For a = M we find an additional feature: The function u'" now has a branch point at a; = muj+. This means that 
there will in principle be a second branch cut in the complex a;-planc. This time, it is the behaviour of the potential 
at r* — oo that leads to the presence of a branch cut. The effective potential near the horizon falls off as 1/r* 
instead of exponentially (as in the non-extreme case). The presence of this second branch cut is discussed further in 
Appendix B. 



III. THE LONG-LIVED KERR QUASINORMAL MODES 



In analogy with standard scattering theory, the quasinormal modes can be interpreted as "resonances" of the black 
hole's gravitational field. This is apparent since they correspond to the poles of the Green's function which, as can 
be seen from (R6t), means that they correspond to 



^out 



= . (41) 



The physical interpretation of this is that the QNMs correspond to waves that are purely outgoing at spatial infinity, 
while at the same time corresponding to purely ingoing waves crossing the event horizon. Given this definition, 
it is straightforward to see that "conservation laws" force the solutions to (^ ) to be complex- valued |Q (in fact, 
for a Kerr black hole, a QNM could exist on the real axis as well, but all attempts to find such modes have failed 
p7| , |3^). It further turns out that there is an infinite set of such mode solutions for each given combination (Z,rn), 
corresponding to complex eigenfrequencies a;„ = w„(M, a, Z,m). The QNM spectra of non-rotating black holes have 
been exhaustively discussed in the literature and are by now well understood In contrast, there have only been a 
few studies of the corresponding Kerr problem. Kerr black hole QNMs were first calculated by Detweiler The 
most detailed calculations were done by Leaver |^ a few years later, and more recently Onozawa ||^ extended Leaver's 
continued fraction approach to discuss near-extreme black holes. Kerr QNMs and excitation coefficients have also been 
calculated by the present authors by integrating suitably chosen phase- functions along paths in the complex -plane 
p6|. These studies lead to the following conclusions: As the black hole spins up, each Schwarzshild QNM (belonging 
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to an infinite family of modes for each I) splits into 21 + 1 separate modes (due to the breakdown of spherical symmetry 
which makes the various m solutions distinct). How the modes are influenced by the black hole's rotation essentially 
depends on whether they are co- or counter-rotating with the black hole. One finds that modes that co-rotate (m > 0) 
with the black hole become longer lived as a increases (i.e. the imaginary part of the QNM frequency is significantly 
reduced). At the same time the oscillation frequency increases. This effect is most pronounced for the I = m modes. 
In fact, as a AI these modes all approach the uj — muj+ point on the real frequency axis, cf. Figure 3 in j^. This 
behaviour can be understood if we note that the / = m modes can be interpreted as being "associated with" the 
black hole's equatorial plane (cf. the symmetry of the spherical harmonics). If such a mode was co-rotating with the 
black hole it would experience maximal frame-dragging effects (see [ p^ for a particularly transparent treatment in the 
eikonal approximation). If we recall that the QNMs can be associated with waves being trapped in the vicinity of the 
unstable photon orbit we see that, as the black hole approaches extreme rotation rates, the co-rotating modes will be 
located closer and closer to the horizon until at a = M when the frequency of these "equatorial" modes is a multiple 
of the rotation frequency of the horizon (a;+). This then makes the extremely slow damping of these QNMs natural. 
Furthermore, it is easy to understand why other QNMs are less effected by rotation. Other co-rotating modes, with 
I ^ m, will be dragged along with the black hole but they never become extremely long-lived (essentially because they 
are not symmetric with respect to the equatorial plane). And the "counter-rotating" modes are much less affected by 
rotation, remaining close to their Schwarzschild counterparts as the black hole spins up. These qualitative results are 
the same for all perturbative fields (gravitational, electromagnetic or scalar). 



A. Approximating the QNM frequencies 

We aim to investigate the extent to which the extremely long-lived QNMs of a rapidly spinning black hole are excited 
by a generic perturbation. As outlined in the Introduction, the main motivation for this study is the question whether 
one should expect these modes to be present in the gravitational- wave signal from astrophysical black holes. If this 
were the case, such signals would be considerably easier to detect than the short bursts that would be characteristic 
of a Schwarzschild black hole. 

As demonstrated by Detweiler the long-lived QNMs can be obtained analytically using an approximation to 
the Teukolsky equation that is valid for a « M and lu « muj-^ . This approximation was first introduced by Tcukolsky 
and Press p^ , and for completeness we reproduce their calculation in Appendix A. Adopting their notation, we define 
new variables 



X 


r ~ 




(42) 


a 


r+ 




(43) 


T 


— M (w — muj- 




(44) 


LU 






(45) 




= - 1/4- 


A , 


(46) 



where X — E + a^uP' — 2maLu. This means that the extreme Kerr limit corresponds to cr — > 0. Furthermore, one can 
show that S is almost purely real (imaginary) for I — m {I ^ m). As can be seen from the equations below, this fact 
distinguishes the I = m case and leads, from a computational point of view, to the presence of the long-lived QNMs. 
From the equations in Appendix A we see that the condition ( |4l| ) corresponds to the following equation; 

T{2iS)T{l + 2iS) 
~ r(-2i(5)r(l - 2i5) 

The left-hand side of this equation has a well-defined limit as a ^ M and lo — > muj^. We represent that limit by 

LHS = ge*x . (48) 

Meanwhile, we see that we cannot have a consistent solution unless t ja —* co as a ^ M. Then the right-hand side 
of the mode-condition can be written (using Stirling's formula) 

RHS = {-%CjTf'^ . (49) 

In other words, a QNM must be a solution to 



r{l/2-2iuj-iS) 



r{l/2 - 2iLd + iS) 



{—2iCja) 



r(l/2 + 2i6j + i5- Air /a) 
T{l/2 + 2ioj-iS-4:iT/a) 



(47) 



f{Lo) = (-8cZ;t)2*'^ - qe'^ = 



(50) 
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Using —Slut = pe"' we see that solutions follow from (remembering that S is real for I = m and u 



p — exp 



c 



X 



2m: 
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From this, we can conclude that there are an infinite number of QNMs such that r — > as a ^ 
black hole we can use uj — ojr^ = ujM, so the long lived QNMs should be well approximated by 



1 

4m 



exp 



X 



2mr 
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This can, of course, be written [B5| 



m 



,(x-2«7r)/25 



4m 



cosC 



,(x-2n7r)/25 



4m 



sinC 



(51) 
(52) 

M. For the extreme 
(53) 

(54) 



It can be verified numerically that sinC > 0, i.e. these QNMs are all damped. 

The approximate result (^4|) provides a useful insight into the nature of the slowly-damped QNMs, but it is only 
relevant for very nearly extreme black holes. The condition ( ^7| ) should, however, remain valid for a larger range of 
spin-rates. Thus we can solve this condition numerically to shed light also on the modes in the near extreme case 
(that may be astrophysically more relevant) . When we do this we find that the / = m QNMs are qualitatively similar 
to ( |5^ ) also for a 7^ M in the sense that, for a give rotation rate, they lie distributed parallel to the imaginary w-axis 
[p6[. As a M the separation between the mode-frequencies decreases until (p3) is recovered. 



B. The excitation of the slowly damped QNMs 



Having established the existence of the long-lived Kerr QNMs we are ready to turn to our main question: Should 
we expect these modes to be excited to a considerable level by a "realistic" astrophysical perturbation? In order to 
begin answering this question, we return to equation (|3^). Let us focus on the contribution from the long lived QNMs 
(and the new branch cut that exists in the a = M case). The remaining contributions to the emerging field should 
not be much altered from the Schwarzschild case. In particular, all the I 7^ m QNMs have frequencies similar to their 
Schwarzschild counterparts and ought to be excited to a comparable level, leading to a signal that dies out after a few 
oscillation periods. The contribution from the "high frequency" arcs can be shown, just like in the Schwarzschild case 

1^, to be zero for late-times, and the field that originates from the familiar lo = branch cut leads to a late-time 
power-law tail which has been recently discussed by Hod For the near extreme case, it follows immediately from 
Cauchy's theorem that the QNM field is 



• +°o r+00 Aout 

$„,(r,0,t)« - ^ / dr:^^e-''^"(*-'---''*)5,,^(0,ac^„)Jim(r:,c^„) 
where we have used the fact that A™{ijj) w (w — w„)a„ near a;„. We have further defined 

Jo 



(55) 



(56) 



Explicit expressions for A™ and deduced for the approximate Teukolsky equation can be found in Appendix A. 
Moreover, the calculation of 



dui 



is straightforward. 

After some manipulations we obtain for the "excitation coefhcient" 



(57) 
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(-_-^-)l/2+2»<i„+»A-„(-2-^ ^4»^„ r(l/2-2iicj„-i(5„) /^27r5„ _ r(l/2 - 2iu!n + i6n)T{l/2 + 2iu}n " i5n) 



r(l/2 + 2iLbn - i6n) 



r(l/2 - 2iu}„ - i6n)T{l/2 + 2itj„ + 16,,) 



xLHS(tJ„) 



duj 



(58) 
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where 5„ = 5(w„). As above LHS(a;) stands for the left-hand side of ( j^7\) but without assuming the double limit 
{u! moj^, a M). The same statement holds for f{uj) — RHS(ijj) — LHS(w). Its derivative at w = a;„ is 



duj 



iRHS(w„) 2r+ -f S'„ 



AM 



^(1/2 + 4ir„/(T„) 



2r+ - 5' - 



AM 



i/)(l/2 + 2iCjn - i5n - AiTn/un) + 2 (5^ ln(-2ic2;„) + 



2iLHS(w„) {(j; (V'(2i,5„) + V(-2i(5„) + ^(1 + 2i5n) + ^(1 - 2i,5„)) 



- {2r++5'^)^{l/2-2iuj„-i5n) + {2r+-5'^)i:{\/2~2iuOn + i5n)} (59) 

where we have used 

dS 



^'-= dc. 



(60) 



and denotes the digamma function [ p6| . 

Given any reasonable initial data we can now "reconstruct" the signal from the long-lived QNMs using (^5|),(|5^) 
and (|59|). Results of such a calculation are presented in Figure |^. 

The equations above are, however, somewhat complicated and it is useful to consider the a — M case for which 
they simplify considerably. Let us first consider the I = m term in (|5^). Bearing in mind that for extreme black 
holes the point uj — muj^ is both a pole and a branch point (see Appendix B), we place the necessary branch cut 
along the family of Green's function poles corresponding to (|54|), see Figure g. After it has crossed the last mode 
(meaning the mode with the largest imaginary part out of the subset of I = m QNMs that approach uj = muj^ the 
cut bends to become parallel to the imaginary w-axis. Choosing an integration contour as illustrated in Figure it is 
straigthforward to show that each of the circles Ci will give a contribution equal to the residue of the pole. That 
is, we get 

I r g-jw(t-r,) / j^out \ 

271- 7branch-cut 2lW V ^ / 



where we have defined 



1 

5m = - y 5™™(0, Afw„)X™™(a;„)e-*""(*-'-*) (62) 



for the mode-sum (with n running from some > to -foo). The quantity £rnm represents the contribution from 
the rest of the branch cut (from point A in Figure |2| downwards). Exphcitly we have 

1 n—ioo 

/:™n = :^e-^"*(*-''-W dpe-'P^'-'-'^lFiuJA + pe^n-FiujA + p)] (63) 
where p = lu — lu^- The function F{uj), defined as, 

= -^Irnm{uj) , (64) 

can be considered as analytic (after analytic continuation) along the integration contour. This will also be true for 
the function F{p) = F{uja + pe'^'") - F{uja + p)- Therefore, we can expand this function in a power series around oja- 

1 p — ioo 

Z:™™ = :^e-'"^ M dp e-'P^'-"-'^ F{Q) + 0{p) . (65) 
Ztt Jo L J 

Integration then yields 

^-uJA{t-rt) 

^mm ^ as t — > +CO . (66) 
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FIG. 2: Schematic description of the QNM location (crosses) and placement of the required branch cuts (solid lines) in the 
case of an extreme Kerr black hole. 



As we will soon argue, this contribution can be neglected at late times. 

We now compute the mode-sum (|62|). The quantity a„ can be deduced from the derivative of the mode-condition 
(or equivalently, by taking the a — M limit in (B3) ). We get 



-SimMSn exp 



1 

26n 



(67) 



Given this, we find that 

^out 



r{l/2 + im~iSn) 



2^s^ T{l/2 - im + iSn)T{l/2 + im — iSn) 



r(l/2 - im - i6n)T{l/2 + im + iSn) 



df 



dio 



Hence, we have shown that 



(68) 



t/<5„ 



(69) 



i.e. each long-lived QNM is excited to an exponentially small amplitude. The total QNM response, ^m, can now be 
approximated by 



Using (p8|) and (Isl), we see that 



'M 



+ 00 



E 



(70) 



'M 



=7V 



^exp -{l^2i5)'^- I3e-^^" 



(71) 



where /? = iuexp(0/2(5 — i(p) / AmM , and we have introduced the retarded time u — t ~ ~ r'^. This can be written as 



Tl = 



5M-^exp -(l-2z<5)^-/3e 



(72) 



with j3 — /Je^'^/*. Since the main contribution to the sum is coming from the large values of n, we can approximate 
the sum by an integral, 



Dm ^ — 

TT 



exp 



-ax — I3e 



dx 



(73) 
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where a = I — 2i5. Expanding the expression in the bracket as a Taylor series and integrating we get 



TT ^— ^ k\(k + a) TT 



after identifying the series representation of the incomplete gamma function 
can deduce that 



Sm 



t 



as t 



oo 



(74) 

Since 7(0;, /?) r(a) as /? ^ 00 we 

(75) 



Since this decays slower than ( p^ ) we can neglect the contribution from the extreme Kerr branch cut. 

Next, we turn to the I > m terms in (|3^). In view of the absence of any poles in the vicinity of a; = mcij+, we simply 
place the branch cut parallel to the imaginary axis. The only multivalued quantity at w = muj^ in the integrand of 
(B§) is the ratio (see Appendix A). A standard branch cut calculation then yields 



$m(r,0,t) 
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— 00 

^out 



{re^'^yM + m/2M) - ——{t/M + m/2M) 



dTe~'''f^' Si„,{9,aT/M + am/2M)Jim{r'^,T /M + m/2M) 



We can split the frequency integral into two parts (for brevity we omit the integrands): 



dT{...} 



dr{...} 



(76) 



(77) 
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The constant e is chosen in order to allow us to use the a — M ^ uj k muj^ approximation 
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2i5 



(78) 



The key point here is that for / 7^ m the quantity 5 is almost purely imaginary near mci;+, with positive imaginary 
part. We then have for r ^ 0; 



2i& 



+00 



(79) 



Hence, the expression inside the bracket in the first integral of the right-hand side of (|77|), will be vanishingly small 
and therefore the contribution from the part of the branch cut close to mw-f is negligible for I ^ m. The second 
integral in (|7^) can be evaluated using arguments similar to the ones involved in the derivation of (|66[). We end up 
with an exponentially decaying field that can be neglected compared to (|75|). 

Finally, we have to consider the contribution from the remaining QNMs. But because all these modes are "short- 
lived" their contribution is important only during the early phase of QNM ringing, cf. the Schwarzschild results 

0- 

The main conclusion of this Section is that even though the contribution of each long lived QNM to the emerging 
field is exponentially small, the total mode contribution cannot be neglected. When summed, the slowly damped 
QNMs of an extreme Kerr black hole give rise to an oscillating signal whose magnitude falls off with time as a power- 
law. This result is particularly interesting since the decay of this signal is considerably slower than the standard 
power-law tail in the Kerr case |l^, ^ |l6j . It is worth mentioning that oscillating power-laws are known to arise in 
standard scattering theory whenever the Green's function has multiple poles js^. In our case, one could argue that it 
is the exponentially small "spacing" between neighbouring poles (see expression (|4|)) that causes the 1/t power-law. 

It is also worth pointing out that the prediction ( |69| ) for the excitation of each individual long-lived QNM agrees 
with the intuitive arguments we made in the Introduction. For example, we see that the leading QNM is not excited 
for a = M (as it corresponds to the limit n — > 00 in (|69[)). This is in agreement with, and generalises, the I — 1 
result (which is valid only for the fundamental mode) obtained in p4| using the eikonal approximation. 
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IV. NUMERICAL RESULTS 
A. Time evolutions 

In the previous Section we used analytic approximations to analyze the long-lived QNMs of a near extreme Kerr 
black hole. We obtained results that agree with our anticipations that these modes would not individually be excited 
to a significant level. But our calculation also provided surprises, the most important being that the large number 
of slowly damped QNMs combine to give a significant signal the amplitude of which falls of as 1/t at late times. 
However, in view of the many approximations involved in the derivation of (|7^) considerable caution is warranted, 
and a confirmation of the analytic prediction is desirable. One way to obtain support for the results would be to 
perform time-evolutions of the Teukolsky equation from given initial data. In other words, the recent effort to develop 
a framework for doing perturbative time-evolutions for Kerr black holes jl^, |ll|, ^ [l^ provides the means for testing 
our analytical predictions. 

We have performed a set of evolutions (for various values of m) using the same scalar field code that was used to 
study superradiance in a dynamical context [T^ . As initial data we have chosen a generic Gaussian pulse originally 
located far away from the black hole. The evolution of this pulse, as it travels towards the black hole and excites 
the QNM ringing is then studied. The numerical results we obtain support the following conclusions: i) In the case 
of extreme Kerr black holes (a = M) we verify the predicted oscillating \/t behaviour for all m 7^ 0, cf. Figure ||. 
ii) For a < M we recover the anticipated exponential fall-off at late times, cf. [|ll|. Still, the emerging signal differs 
considerably from a single QNM oscillation at intermediate times for near extreme black holes. This agrees well with 
the analytical results from the previous Section: For near extreme black holes a large number of long-lived QNMs are 
excited to roughly the same level, and for a considerable time window the resulting signal is a superposition of many 
slowly decaying exponentials (each with a very small amplitude). As we depart further from a ~ M we retain the 
standard result: The signal is completely dominated by the slowest damped QNM. iii) For axisymmetric perturbations 
(to — 0) the numerical evolution recovers the standard power-law tail. For our particular choice of initial data (that 
contains the / = multipolc) the tail falls off as . This agrees with the predictions of, for example, Ori and Barack 
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FIG. 3: A numerical evolution showing the late-time behaviour of a scalar field in the geometry of a rapidly rotating Kerr black 
hole. We show (on a logarithmic scale) the field as viewed by an observer situated well away from the black hole for a = M. 
At late times the field falls ofi^ according to an oscillating power- law with the amplitude decaying as 1/t. The data corresponds 
to a narrow Gaussian (initially centered at r* = 50M) that hits the black hole, and is observed at r, = IQM 

Although our numerical simulations generally support the analytic results, there is still room for some caution. It is 
very difficult to investigate the late-time behaviour of a perturbed Kerr black hole using the Teukolsky code. After all, 
we would like to be able to distinguish a slowly damped exponential from an oscillating power-law at very late times. 
Given our current numerical code we cannot obtained absolute proof of our analytic results. Detailed convergence 
tests show that the code has the expected properties, but also unveil that it is difficult to reliably determine the late 
time oscillating tail predicted by our analytical work. While the overall amplitude envelope can be determined for 
many hundred dynamical timescales as the resolution is increased, the signal typically goes out of phase on a shorter 
timescale. This means that, even though the numerical results provide support for our theoretical predictions, one 
would have to build an evolution code that could be trusted at very late times, eg. based on a double-null evolution, 
in order to obtain definite results. 
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FIG. 4: We compare the late-time QNM oscillations in a typical scalar field evolution (initial data and observer location are the 
same as in Figure ^ for m = 2 and a = 0.99Af to the damping rates of the two slowest damped QNMs (indicated by dashed 
lines). From this figure it is clear that the late-time behaviour is not well represented by either of these modes individually. 

B. Reconstructing the long-lived QNM signal 

In addition to investigating the late-time field through direct numerical evolutions, we can reconstruct it from our 
analytical results. Using expressions ( |55| ) - (^9|) it is a straightforward task to calculate the long-lived QNM signal 
provided that we have already acquired the corresponding mode-frequencies. In this way one can expect to obtain a 
fairly accurate description of the full signal for sufficiently late times (i.e. after some e-folding times of the dominant 
short-lived QNM: typically the slowest damped m = —I mode), spanning a time window of several hundred dynamical 
timescales. Eventually one would expect this signal (after all, it decays exponentially) to give way to the standard 
power-law tail. However, as is clear from all our illustrations this will not happen until at very very late times. 

In Table | we present a small subset of our numerical data for long-lived QNMs frequencies and the corresponding 
excitation coefficents for near extreme black holes. As anticipated, these amplitudes vanish as a — s- M. Moreover, 
as already pointed out, the higher overtones are excited to a level comparable to the leading mode. In Figure |^ we 
illustrate the typical field as constructed from the long-lived QNMs. For this particular calculation, we have chosen, as 
initial data, a Gaussian pulse with angular dependence ^ 122(6*, </>). It follows then from (p5|), that only the I = m = 2 
field component will be important at late times. The result in Figure ^ should be compared to the fully numerical 
evolution data in Figures |[^. 

V. A PHYSICAL INTERPRETATION: THE SUPERRADIANCE RESONANCE CAVITY 

Our analytical results, supported by the outcome of numerical evolutions, give rise to an intriguing picture. It would 
seem as if the extremely long-lived QNMs will, even though they are hardly excited at all individually, dominate the 
signal from a very rapidly spinning black hole. The emergence of this late-time signal, to which many QNMs contribute 
would be a new phenomenon in black- hole physics. Given that both the standard QNMs and the power-law tail have 
simple intuitive explanations it may be worthwhile trying to understand the extreme Kerr case in a similar way. From 
the evidence provided by our numerical evolutions, see Figure ^ we propose the following "explanation" . Consider 
the fate of an essentially monochromatic wave that falls onto the black hole, cf. Figure || for a schematic description. 
Provided that the frequency is in the interval < w < muj^ the wave will be superradiant: scattering of these waves 
in the black hole's ergosphere results in their amplification by extraction of the black hole's rotational energy. In 
effect, this means that a distant observer will see waves "emerging from the horizon", cf. (^2|), even though a local 
observer sees the waves crossing the event horizon (at r+) |2^. In addition to this, one can establish that the effective 
potential has a peak outside the black hole (which is not immediately obvious since the "potential" Q is frequency 
dependent in the Kerr case) for a range of frequencies including the superradiant interval. Indeed, we can find the 
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FIG. 5: The signal (shown on a logarithmic scale) from the long-lived QNMs (the first ten modes are included here), for an 
a = 0.9999M Kerr black hole. The required mode-data is given in Table pi The observer is located at r» = lOOM, 9 — 11/2 and 
a narrow Gaussian pulse (centered at r, = 50Af) was used as initial data. The dashed line represents the damping rate of the 
slowest decaying QNM. It is clear that for several hundreds of black hole's dynamical times, the slowest damped mode is not 
dominant. This is in contrast to the results for slowly rotating black holes. 



following approximate expression for the potential near the horizon (valid for a — M only) 

Q(r, Lo)^iLO- mu+)^ - 2(c^ - mio+)- + O ( ^Jl^i^lIMl) (80) 

r* \ ri J 

This shows that, for uj « muj^^ there will be a peak (corresponding to a minimum of Q) just outside the horizon. 
This can be verified by graphing the exact potential (|2^) , cf . Figure |^. 




FIG. 6: Schematic explanation of the new phenomenon seen in the numerical evolutions of Kerr perturbations. The left panel 
illustrates the standard scenario: An infalling pulse excites the QNMs that then propagate to infinity and the horizon. At late 
times, backscattering due to the curvature in the far-zone dominates and leads to the familiar power-law tail behaviour. Right 
panel: For frequencies that i) lie in the superradiant regime, and ii) experience a "potential peak" in the region [r+,cx3] there 
will be a superradiance resonance cavity outside the black hole. At late times, the waves leaking out of this cavity to infinity 
dominate the signal. 



The combination of the causal boundary condition at the horizon effectively corresponding to waves "coming out of 
the black hole" (according to a distant observer) and the presence of a potential peak leads to waves potentially being 
trapped in the region close to the horizon. In effect, there is a "superradiance resonance cavity" outside the black 
hole. Again according to a distant observer, waves can only escape from this cavity by leakage through the potential 
barrier to infinity. Since the superradiant amplification is strongest for frequencies close to rna;+, waves in the cavity 
experience a kind of parametric amplification and at very late times the dominant oscillation frequency ought to be 
mw+. This is exactly what we have deduced from our analytic and numerical calculations. In Figure ^ we present 
a series of snapshots of the evolution of a scalar field around an extreme Kerr black hole. In this series of pictures 
one can see how "trapped" oscillations develop in the vicinity of the black hole's horizon. After an initial drop in 
amplitude (roughly at the timescale of the e-folding time of the rapidly damped "Schwarzschild-like" QNMs) these 
oscillations decay very slowly. These snapshots indicate the presence of a standing wave in the region just outside 
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the horizon. In the extreme black hole case leakage from the superradiance cavity leads to the observed 1/t decay. In 
the near extreme case, the existence of the cavity provides an intuitive explanation for the extremely slow damping 
of corotating QNMs with frequencies close to maj_(-. In a way, the superradiance cavity we have described here can 
be viewed as the black hole analogue of the potential "well" present in the gravitational field of ultracompact stars 
p8[ which also leads to the existence of a long-lived family of w-modes p9] . 
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FIG. 7: Graph of the I — m = 2 effective radial potential Q{r,,uj) as a function of r*, for lu = 0.9mLj+ (dashed curve), mu!+ 
(solid curve) , l.lma;+ (dotted curve). The black hole spin is a = 0.999M. Note the appearance of a potential barrier {Q < 0) 
as soon as we enter the superradiant frequency regime uj < mui+. 



VI. CONCLUDING DISCUSSION 



We have presented the results of an investigation into the late-time behaviour of a perturbed Kerr black hole. An 
analytic calculation for scalar fields in the geometry of an extreme Kerr black hole provided two important results. 
The first concerns the level of excitation of the QNMs that become very slowly damped as a ^ M. We find that 
these modes are much more difficult to excite than their rapidly damped Schwarzschild counterparts. This means 
that (individually) these modes may not be easy to detect with the new generation of gravitational- wave detectors. 
However, there may still be a detectable, slowly damped, QNM signal owing to the fact that a large number of 
virtually undamped QNMs exist for each value of m 7^ 0. We find that these modes combine in such a way that the 
field oscillates with an amplitude that decays as 1/t at late times. This decay is considerably slower than the standard 
power-law tail. We have used numerical time-evolutions of the Teukolsky equation to verify this analytic prediction 
for extreme black holes. 

After extending our study to near extreme black holes we find that, even though the QNM signal decays expo- 
nentially in the familiar way, there is still a considerable time interval where a large number of small-amplitude 
slowly-damped QNMs are present in the signal. This means that the signal from a rapidly spinning black hole would 
not be well represented by a single QNM approximation. Only at very late times do the leading mode become 
dominant and eventually it should give way to the power-law tail predicted in previous work p^ . 

Whether our results are astrophysically relevant or not is in many ways still an open question. In particular, the 
excitation of the long-lived QNMs by more realistic initial data (representing, for example, the close-limit approxima- 
tion of merging Kerr black holes) must be studied. One should obviously also try to understand whether one would 
expect to find almost extreme black holes in Nature (recall that Thorne has shown that accretion cannot spin a 
black hole up beyond a = 0.998M). It seems quite reasonable, however, to expect that the effects we have observed 
will play a role at intermediate times for rapidly rotating non-extreme black holes. If our predictions are correct one 
must investigate in detail to what extent the late-time signals from a rapidly rotating black-hole are detectable even 
though they contain a large number of QNMs (each with a small amplitude). In particular, it is crucial to determine 
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FIG. 8: A series of snapshots of a typical numerical evolution of a rapidly spinning black hole. We show the scalar field a 
function of both r, and 9 at various times (as indicated in the upper right hand corner of each image). The images illustrate 
how "trapped" oscillations develop in the vicinity of the black hole's horizon. After an initial drop in amplitude, cf. the first 
four images, (roughly at the timescale of the e-folding time of the rapidly damped "Schwarzschild-like" QNMs) these oscillations 
decay very slowly. These last four snapshots indicate the presence of a very slowly damped standing wave in the region just 
outside the horizon. 



whether one can hope to infer the black-hole parameters from a signal that contains such a superposition of QNMs 
(perhaps using techniques similar to those introduced in |2^). It is also interesting to ask whether there exists a 
critical value of the rotation parameter a above which the new effect we have observed becomes relevant (recall that 
our approximate modes are only relevant for a « M) . More detailed numerical work is needed to answer this question. 
Finally, we note that the observed phenomenon can be intuitively explained in terms of a "superradiance resonance 
cavity" outside the black hole, in which waves of certain frequencies are effectively trapped. As the waves slowly leak 
out from this cavity, they give rise the "long-lived" QNM signal observed at infinity. In conclusion, it is interesting 
to note that, even though black-hole perturbations is a very well researched field, it can still provide interesting and 
surprising results. 
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APPENDIX A: APPROXIMATE SOLUTIONS OF THE TEUKOLSKY EQUATION FOR NEAR 

EXTREME ROTATION 

In this Appendix we reproduce the solution of the Teukolsky equation for a ~ M, uj « mLO+, as given by Teukolsky 
and Press Throughout this paper we have (as far as possible) used the notation of their Appendix Al. The 

equation to be solved is not the "usual" Teukolsky equation, but the one that follows if we consider the problem in 
"ingoing Kerr" coordinates. These are related to the Boyer-Lindquist coordinates by 

v^t + j—^dr, (Al) 

^ ^ ip + J ^dr . (A2) 

(A3) 

In these coordinates the separation of variables follows from (note that the radial functions here differ from the one 
used in Section IIA by the factor (r^ + a^)^^/^) 



^im = j Rim{u;,r)SUO,aLu)e'"'^-''''^du; . (A4) 
Comparing this expression to the standard Boyer-Lindquist decomposition 

<^im = J RimiuJ, r)SUO, auj)^"'^''^'dw , (A5) 
we can immediately relate the radial wavefunctions via 

-Rzm(w, r) = Rim{(jJ, r) exp 

Since 



(A6) 



K 



A 1 fc?"* , as ^ — oo . ^ ' 



where we recall that k = lo — muj+. Hence, we see that the "physically acceptable" solution behaves as 

1 as r ^ r-|_ 

e^*"''* + Z'"r ^ as r 



In ingoing Kerr coordinates, the radial Teukolsky equation takes the form, 

x{x + a)^-^ - {2iu;x^ + 2x{2il^ - 1) - (t + 4ir}^^ - {2iw(x + 1) + A}i?/,„ = (A9) 
The variables x etcetera were defined by (E^)-(E^). The double limit a — > Af, lo — > mujj^ corresponds to cr 0, 



T — > 0. Let us first consider this equation in the limit when x » max (cr, r), i.e., for large radii. Then (A9) is well 
approximated by 



^'fRlrn _ ^2-^^2 _^ 2x{2iuj - 1)}^^ _ {2tu,{x + 1) + A}i?i™ = . (AlO) 



da; 2 dx 
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A solution to (AlO) that satisfies (A8) can be written in terms of confluent hypergeometric functions 

Rl„^ = Ax-^/^+^'^+'^M{l/2 + 2iLj + i5,l + 2iS, 2iCox) + B{5 -5) 



(All) 



where A,B are constants and the notation [6 —5) means "replace 5 by —5 in the preceding term". We next turn to 
the case when x << 1, i.e., try to find a solution that is valid close to the black hole's horizon. The original equation 
simplifies (slightly) to 



x{x + cj)'^^ - {2x{2iLo + Air)}^ - {2iCo + A}i?,,„ = . 



dx'^ ^ ■' " dx 

This is the hypergeometric equation, and one solution can be written as 

Rhn = 2i^i(l/2 - 2iuj + i(5, 1/2 - 2iuj -i5,l- Air /a, -x/a) 



(A12) 



(A13) 



It is straightforward to verify that 2-^1 ^ 1 as x — > 0, which means that this solution has the desired "purely ingoing 
wave" behaviour close to the event hori zon. The solutions ( All ) and ( A13 ) can be matched in the overlap region, 
max((T,T) < x < 1. The x ^ hmit of ([All|) yields, 



im 



Ax-^/'^+'^'^+'^ + B{5 ^ -5) 



Similarly, for x ^ 00, ( A13 ) becomes, 

r(l - AiT/a)T{2i5)a^/^-'^''^-''^ 



Rlr. 



r(l/2 - 2iuj + i5)T{l/2 + 2iCj + i5 + Air/a) 
We can extract A and B by matching the solutions ( |A14| ) and ( |A15| ), 

r(l - 4ir/(T)r(2i5)crV2-2»<i~i5 



A = 

B = 



r(l/2 - 2iCj + i5)T{l/2 + 2iuj + i5 - Air jo) 
r(l - 4ir/(T)r(-2i(5)(7i/2-2»c:.+i5 

r(l/2 - 2ia) - i(5)r(l/2 + 2iLj - i6 ~ Air/a) 



On the other hand, approximating ( All ) for x ^ 00 we get for the amplitudes 

Z 

r{i + 2iS) 



1(1/2 — 2iuj + 10) 



r+A 



T{l/2 + 2iu) + i5) 



— (2^a>)-l/2+2^.i-^<5 ^ ^ 



(A14) 
(A15) 

(A16) 
(A17) 

(A18) 
(A19) 



Using (|A6| ), it is easy to see that A'" = and A°^* = where and A™ are the asymptotic amphtudcs in 
the standard Boyer-Lindquist decomposition. Explicitly we get 



A'^'ioj) ^ r+(-2it2))-i/2-2«^-«5^i/2-2»^-«5 



T{-2iS)T{l - 2id)Til - Air/a) 
(r(l/2 - 2iL2; - i5))^T{l/2 + 2iuj + iS - Air/cr) 



r{2iS) T{l + 2i6) /r(l/2-2icj-i5) 
T{-2iS) r(l - 2i5) \r{l/2-2i^ + iS) 



+ {-2iLua) 



r(l/2 + 2iu} + iS- Air/ a) 



A°"*(w) = r+{2iLb) 



V{l/2 + 2iLj-i5- Air /(j) 

T{-2i5)T{l-2i5) T{l~AiT/a) 
r(l/2 + 2iuj - i(5)r(l/2 - 2iCj - i5) r(l/2 + 2iuj + iS- Air/ a) 



(A20) 



r{2iS)T{l + 2iS)T{l/2 + 2iu} - iS)r{l/2 - 2iLo - iS) ^2^^ ^(1/2 + 2iL2) + iS ~ 4iT/cr)l 

r(-2i(5)r(l - 2id)T{l/2 - 2iu + i(5)r(l/2 + 2iuj + id) ^ ^ T{l/2 + 2iw - id - Air/a)} ' 



and their ratio will be 

^out 



_,)l/2+2.^+..(2^^)4.^r(l/2-2zC^-Z^) 

' ^ ' T{l/2 + 2iuj-i5) 

T{2i5)T{l + 2i5)T{l/2 + 2iuj - i5)T{l/2 - 2iuj - iS) 
r(-2i(5)r(l - 2iS)r{l/2 - 2iuj + i5)T{l/2 + 2iu) + i5) 



, T{ll2 + 2iQ^ + z5-AiTla) 

+ [ZlUJa) — : ; ; — - 

^ ' T{l/2 + 2iCj-i5-AiT/a) 



T{2i5) T{l + 2i5) (T{l/2-2iLo -i5) 



T{~2i5) r(l - 2i5) \T{l/2 - 2iuj + id) 



{—2iLba) 



r(l/2 + 2iu} + id- Air/ a) 
T{l/2 + 2iiv-i6- Air/a) 



(A22) 
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In the special case of an extreme Kerr black hole we have ct = 0, and ( A22| ) then reduces to 

^out 
in 



A" 



T{2i5)T{l + 2i5)T{l/2 + 2iujM ~ i5)T{l/2 - 2iujM - id) 
r(-2i(5)r(l - 2i6)T{l/2 - 2iujM + iS)r{l/2 + 2iLuM + iS) ^ 

2 

,2iS 



T{2id) T{l + 2i6) f Til/2 - 2iuiM -i5) 
T{-2i5)T{l-2i5) \ T{l/2-2iujM + i5) 



-1 



2iS 



(A23) 



Note that this expression is multivalued at w = mw+ due to the presence of the t term (see Appendix B 



APPENDIX B: ANALYTICAL PROPERTIES OF SOLUTIONS OF THE TEUKOLSKY EQUATION 

In this Appendix we study the analytical properties of the u'" function in the complex w-plane. A much more 
rigorous and detailed treatment has been provided by Hartle and Wilkins However, as we will show here their 
conclusion that u™ has a branch point for a < M is incorrect. Instead, we find that u™ has a series of (physically 
insignificant) poles. Only for the special case a = M do these poles disappear, and u™ has indeed a branch point. 

To begin we write the radial Tcukolsky equation (|2^ ) in a more compact form (in this Appendix we shall be 
considering perturbations of an arbitrary spin s field) 

^ +Q{r,,iu)u = . (Bl) 

It is a well-known fact in scattering theory that certain analytical properties (with relevance for the late-time 
behaviour of the field) of the solutions to an equation of the form (Bl), can be found by studying the asymptotic 
behaviour of these solutions as r* ±00. Therefore, we consider (Bl) in the limit r r+. For a < M the potential 
has the asymptotic form, 

(^(r,,^) « (cj - mw+)^ - T4(^) + V"o(w)e'""* as — > -cx) , (B2) 

where 

and c = (r_|_ — r_)/2Ai^r_|_. The explicit form of the Vo(lo) function will not be required in the following. Setting 



equation (Bl) becomes 

Ke-] u = , (B5) 

near the horizon. This equation can be solved by employing the standard Born approximation. However, it turns out 
that we can solve it exactly |]37| . By setting y = e'^''* we can rewrite the equation as 

y T^ + y^ + -i^ -'^oy)u = o (B6) 

cty ay c 



This is a Bessel-type equation and the solution u'" is 

/ 1/2 \ ^'^/'^ / 1/2 

«'"(r.,a;) = r(l - 2^^l/ c)J_,,n,c ^^^^e^'/^ ) (B7) 



From the small argument approximation for the Bessel function (recall that r» ^ — cx) as we approach the horizon), it is 
easy to see that this solution indeed satisfies the appropriate boundary behaviour at the horizon: u'" — > A~''/^e~*'^''* . 
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The presence of the F function signals the existence of simple poles at the points 1 — 2iVl/c = — n (where n a 
non-negative integer) or, equivalently, at frequencies 

w = + i(s - n- 1)^^^^— — (B8) 

If we had chosen to use the Born approximation, then at leading order, we would had picked up the first of these 



poles. The existence of these poles can be directly deduced by inspection of equation ( A13 ). These poles have, 
however, no physical significance. The black hole Green's function does not inherit them since the factor r(l — 2iVl/c) 
will be cancelled by an identical fa ctor in the W ronskian (in the r» +oo form of u™ these poles are contained in 



the amplitudes A™, A™*, see eqs. (A16), (A17) of Appendix A ). As a consequence, these are called "false" poles in 
quantum scattering theory fs?" 



Due to the appearance of the Bessel function of a noninteger power in (B7), one would expect m'" to have branch 
points at those frequencies Wb that solve Vo{uJb) = 0. However, since J^{ze'^^ = e'^'^'^^ J^{z) Q we are effectively left 
with a single-valued function. It is a quite general result in scattering theory that exponentially decaying potentials 
cannot lead to the existence of branch points |37) . 

We now turn to the a = M case. Approximating the effective radial potential near the horizon for a ~ M yields 

Q{r,iLj) fa {lu ~ TOa;_|_)^ — 2{lj — mLJ-^-) as — s- — oo . (B9) 



Near the horizon then, eqn. (Bl) effectively becomes 



||+(ff + |)..0 (BIO, 

with Cl ~ u! — mu!^ and (3 — — 2(to — is)tt. This equation can be easily transformed to a Whittaker equation 
with two independent solutions WK.p,{2iflr^,), W_k,/j(— 2if2r,) where k = im + s, ii = ±1/2. The Whittaker function 
VF„,^(z) is defined by, 

W^.^,{z) = zi/2+'^e-^/2[/(i/2 + /i-K,l + 2^,z) , (Bll) 

(B12) 

Taking a large argument approximation for the confluent hypergeometric functions we deduce that WK^fj.(2iQr^) is the 
solution with the desired "ingoing wave" behaviour at the horizon. Clearly, WK,fj,{z) is a multi- valued function with 
a branch point at z = 0, which in the present case corresponds to w = ma;+. This is in agreement with the result of 
Hartle and Wilkins jss) . In this case the branch poi nt is inherited by the relevant Green's function. More specifically. 



it is contained in the ratio see equation (A23) 



APPENDIX C: GRAVITATIONAL PERTURBATIONS 

So far, we have mainly considered the model problem of a massless scalar field. We did this for reasons of simplicity. 
However, it turns out that our results are easily extended to the physically more interesting case of gravitational 
perturbations. In fact, the calculation proceeds very much along the same lines. This means that all conclusions 
reached in the main part of this paper remain valid also for the gravitational case. 

In the Teukolsky formalism, gravitational perturbations are represented by the the scalar functions with s = ±2 
p2| . As before, we can separate the dependence on the (p angle 

= A-V2(^2^^2)-i/2 y ^^^{r,e,t)e'^^ (CI) 



m—-~oo 



Further decomposition is possible in the Fourier domain 

p-\-oo-[-ic 



, I _|N J — OC + ic 



Z— max( |m| , |s| ) ' 

where sRhn solves (pTi) with the potential. 



, , -2is{r~ M)K +UiLurs-2amuj-a'^uj'^ - E)A dG ^ , , 

^("'") = FT^^p rf^-G (C3) 
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with G = s{r — M)/(r'^ + a^) + rA/(r^ + a^)^. The sohitions of the homogeneous equation ( p7| ) which are used for 
the construction of the Green's function are now 



and 



^-Sgjojr, a,s r ^ +00 , 



A */^e as r ^ r+ , 



Repeating the same steps as in Section IIA, we end up with an expression similar to ( 38 ) 

^out 



/Trr^ I /II,) 



/— inax(|m|,|s|) 



— 00+ic 



(C4) 



(C5) 



(C6) 



The analytical properties of the Green's function are identical to the scalar case (recall that the analysis in Appendix 
B is valid for an arbitrary spin field). The long-lived QNMs are now solutions of 

V{2i5)T{l + 2i5)T{l/2 + s~2iCj-i5)T{l/2-s-2iu)-i5) _ ^^^T {I / 2 ^ 2iuj + i5 - Air / a) 



T{-2i6)T{l~2iS)r{l/2 + s~2ii:u + iS)r{l/2-s-2iLu + i6) ' ' T{l/2 + 200 - iS - Air/a) 

This equation turns out to be invariant under the change s — )■ — s. It thus follows that the s — ±2 components share 
the same set of long-lived QNMs. 

We need to find approximate (a « M, lo « ma;+) expressions for the coefficients ^A™, ^A""'. Press and Teukolsky |2^ 
give the relevant approximate radial wavefunctions in a Kerr-ingoing coordinates. We need to relate their solutions to 
the respective wavefunctions s^'", su"^ in order to extract the asymptotic amplitudes and In the standard 

Boyer-Lindquist coordinates, a spin-s field can be decomposed as 

= ,7^^,„(r,c^),5^r(^)e*'"^-^"* . (C8) 

Similarly, in Kerr-ingoing coordinates, we get 

= ,7^^^(r,c^),5^r(^)e™^e-'"'^ . (C9) 

These two representations are related by 

(y) . (CIO) 

We adopt the following normalisation for the "in" radial wavefunctions, at r — > 00 (hereafter, for brevity, we drop the 
Im subscript) 

T-)in r A ''e as r — ^ r-|_ , (Tnl 



1 as r — + r+ , 



as r ^ -l-oo 



(C12) 



We are mainly interested in the s = — 2 component, as it leads to the dominant contribution to the gravitational 
waves reaching infinity Inspection of (C5) and ( |C11 ) gives. 



-2B' 



(C13) 



From (CIC) and (Oil) we obtain. 



/A227^'" exp 



—dr 
A 



(C14) 
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where / some proportionality factor. Furthermore, it is possible to relate the two s = —2 fields, by means of the 
Starobinsky identity |27|: 



with X = E + a^Lo^ — 2aujm. Combining (C14) and (C15) we get 



-27^ = / 



A2 d 



4|C|2 dr 



(C15) 



where J+ = d/dr + iK/A and 

|Cp = (A^ + iaium - 4:a^Lj^)[{\ - 2f + 36awTO - 360^^^] + (2A - \)[<d'6a^u? - 48awTO) + \^^u?{M'^ - a^) (C16) 



Eqns. (IcTil ) and (|cT7| ) can be used to relate the coefficients (-aA''^, -2^°"*) to 2^°"*) and (_2^" , -z. 

respectively. The first of these connecting formulas simply yields _2^"^/-2^°"* = 2-^'"/2^°"*- The QNM frequ 
therefore satisfy 



(C17) 
lencies 



2A" 



^out 



(C18) 



This equation yields (C7) for s = 2, which is, as pointed out, invariant under the operation s — > — s. It is easy to 
verify that the long-lived gravitational QNMs are given by an expression similar to (p4). We have. 



-(^) = 



r^^(2ia;) 



r(-3/2-2iLj-(5) 
r(5/2 + 2iQj - iS) 
T{2i5) r(l + 2i(5)r(5/2-2ia;- 



i5)T{-3/2-2iLd~iS) 



T{-2i5) r(l - 2i5) r(5/2 - 2^0) + i5)T{-3/2 - 2iuj + iS) 
T{2i5)T{l + 2i(5)r(5/2 + 2iuj - i5)T{^/2 - 2iLj - id) 



+ {~2iuja) 



2iS 



T{l/2 + 2ii:u + iS~4iT/a) 



T{-2iS)T{l - 2i5)T{5/2 - 2iu! + iS)r{5/2 + 2i(b + iS) 



{2iLua) 



r{l/2 + 2iu} - iS - Air /a) 
r(l/2 + 2iLb + i6- -iiT/a 



r(l/2 + 2iu} -iS- Air/a) 



(C19) 



It follows that the analysis presented in Section HE can be used also in the gravitational case, yielding results indentical 
to the scalar ones. 

We note that the grav itational perturbation problem can alternatively be approached exclusively via the s — —2 
fields, see equation (C17). This choice was made by Sasaki and Nakamura |^ in their study of the late-time response 
of an extreme Kerr hole perturbed by an orbiting test particle. Their results do not agree with ours in that they find 
that the late-time signal decays as ^ Such a power-law would lead to an (unphysical) divergent flux integral 

p5| . We believe that the calculation of Sasaki and Nakamura contains an error in the relation between the relevant 
asymptotic amplitudes. According to their result, the ratio of the asymptotic amplitudes _2-S"V-2-S°"' introduces 
additional powers of r in the integrand of (|C6| ) which then leads to the different power-law behaviour. We have 
perfomed the calculation (using equations ( CI?] ), (All)) and found 



(C20) 



with 



D{uj) = {2i6jr' ( {26, + Sf + ^{2u; + + ^ j ^2^2, - S)^ + ^(2c2; -3)^ + ^ 



(C21) 



which differs crucially from the result of Sasaki and Nakamura. However, given that several derivatives of the 
asymptotic solutions are required to derive this result, and that significant cancellations occur, it would be very easy 
to make a mistake in this calculation. 
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n Rc cj„M Im u)„M Rc (MA°"Vq„) Im (AfA°"Va„) 






0, 


.99324 


-0, 


.00341 


0.00216 


-0.00369 


1 


0. 


.99322 


-0, 


.01020 


-0.00691 


-0.00928 


2 


0, 


.99321 


-0, 


.01699 


-0.02167 


0.01501 


3 


0. 


.99320 


-0, 


.02385 


-0.02167 


0.01501 


4 


0, 


.99317 


-0 


.03067 


-0.01288 


0.03089 


5 


0, 


.99313 


-0 


.03749 


0.00398 


0.04030 


6 


0, 


.99308 


-0 


.04432 


0.02439 


0.040882 


7 


0, 


.99303 


-0, 


.05115 


0.04427 


0.03253 


8 


0. 


.99298 


-0, 


.05798 


0.06043 


0.01643 


9 


0. 


.99291 


-0, 


.06482 


0.07055 


-0.00550 



TABLE I: Frequencies and excitation coefficients for the first ten I — m — 2 QNMs, for an a = 0.9999M Kerr black hole. These 
results illustrate the fact that many QNMs are excited to a comparable amplitude (cf. the excitation coefficients listed in the 
last two columns). 
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